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ABSTRACT 

We investigate the formation and evolution of giant molecular clouds (GMCs) by 
the collision of convergent warm neutral medium (WNM) streams in the interstellar 
medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on 
the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the 
mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds 
formed by supercritical inflow streams proceed directly to collapse, while clouds formed 
by subcritical streams first contract and then re-expand, oscillating on the scale of tens 
of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 
4 /iG show that only supercritical or marginal critical streams lead to reasonable star 
forming rates. This result is not altered by the inclusion of ambipolar diffusion. 3) The 
GMCs M2FR is a generally increasing function of time, whose growth rate depends on 
the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly 
fiuctuating function of position in the clouds. This implies that a significant fraction of 
a cloud's mass may remain magnetically supported, while SF occurs in the supercritical 
regions that are not supported. 5) In our simulations, the SFE approaches stationarity, 
because mass is added to the GMC at a similar rate at which it converts mass to stars. 
In such an approximately stationary regime, we find that the SFE provides a proxy of 
the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy 
of the low-M2FR regions within the gravitationally-contracting GMCs, so that the 
latter naturally segregate into a high-density, high-M2FR "core" and a low-density, 
low-M2FR "envelope", without the intervention of AD. 
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1 INTRODUCTION 

Magnetic fields in molecular clouds (MCs), and the gradual 
redistribution of magnetic flux within them by ambipolar 
diffusion (AD), have been thought to be crucial ingredients 
in regulating star formation (SF) and its efii ciency (SFE) for 
over two decades (see, e.g., the reviews bv lShu et al1ll987l : 
[Mouschoviaa. 1 99 L . and references therein). In early studies, 
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MCs were considered to generally have strongly magnet- 
ically subcritical mass-to-flux ratios (M2FRs ), and super- 
critic al clouds were believed to be rare (e.g., iMouschoviad 
Il99ll . sec. 2.1) so that the time for AD to allow their 
cores to become supercritical turned out to be very long, 
of order 10-20 times larger than the clouds' free-fall time 
(tff) (e.g., ICiolek fc Mouschoviaslll994l : lBasu fc Mouschoviad 
Il994h . In this "standard model" of magnetically-supported, 
AD-mediated subcritical clouds, the general mechanism for 
the formation of low-mass stars was through the slow gravi- 
tational contraction of isolated cores containing a very small 
fraction of the clouds' mass, thus accounting for the very low 
observed global SFE o f giant MCs (GMCs) (|Mvers et all 
1 19861 : lEvans et al.ll2009l ) . 
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H owever, subse quent studies have s ugRested that both 
MCs llMcKed ITosgl') and their clumps (ll3ertoldi fc McKeel 
1992"; Mvcrs & Goodman'l988VCrutcher"l999';'Bourkc ct al.' 
2001; Crutchcr ct al. 2003; Troland fc Cru tchcr. .2008 ) arc 
close to being magnetically critical, with a moderate pref- 
erence for being supercritical. This implies that, if a MC is 
subcritical, it is expected to be only moderately so as well, 
in which case the time for cores within it to become locally 
super critical may be almo s t as short as the cores' free-fa ll 
time l|Ciolek fc Basull200ll : IVazauez-Semadeni et al.ll2005l ). 
Moreover, MCs are generally beli eved to be supersonic ally 
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and references therein), and in this case, AD must be treated 
nonlinearly. As a consequence, its characteristic timescale 
is expected to decrease, suggesti ng again that it rnay be 
comparable to the free-fall time (jPatuzzo fc Adain^ l2002l : 
iHeitsch et al.ll20o3 ). Thus, the gravitational contraction and 
collapse of a star-forming region must occur rapidly, essen- 
tially in the timescales corresponding to those of a mag- 
netically supercrit ical region, which i s of the order of a 
few free-faU times JOstriker et ai1ll999l: iHeitsch et al.ll200ll : 



IVazauez-Semadeni et al.ll2005l : ICalvan-Madrid et al.H2007h . 
although the subcritical envelope may still be held up by the 
magnetic field, since, at the envelopes' lower typical densi- 
ties, the AD timescale is indeed much longer than the dy- 
namical timescales. 

Moreover, the realization that the majority of MCs may 
be supercritical, and that most stars , including low-mas s 
ones, form in cluster-forming regions (|Lada fc Lada ''2003'). 
have forced a reconsideration of the problem, to accomo- 
date the fact that the standard model's paradigm of low- 
mass star formation in strongly subcritical clouds may be 
the exception rather than the rule, even for the formation of 
low-mass starsQ Unfortunately, if mo st MCs are supercriti- 
cal, on e must once again face the old IZuckerman fc Palmeil 
l|l974h conundrum that the Galactic star formation rate 
(SFR) should be roughly two orders of magnitude larger 
than the one obser ved, of ~ 3-4Mp) yr^^ (see the supple- 
mentary material of lDiehl et al.|[200a ). This is because glob- 
ally supercritical clouds cannot be supported by the mag- 
netic field, and thus should be collapsing as a whole. 

Turbulence is often invoked as an additional source 
of support against the clouds' self -gravity, as if it wer e 
simply an extra source of pressure l|Chandrasekhadll95ll ). 
Such a treatment, however, neglects the fundamental prop- 
erty of turbulence that the largest velocities occur at the 
larges t scales, a property which is refiected in iKolmogorovl 
l|l941l )'s famous energy spectrum of the turbulence. Stud- 
ies taking this into account have generally only considered 
it from the point of v i ew of the energetics involved (e.g . , 
iBonazzola et al.l Il987l : IVazguez-Semadeni fc Gazoll Il995l ). 
but have neglected the vector nature of the velocity field. 
Having the largest velocity differences at the largest scales 
implies that the effect of classical vortical turbulence within 
a cloud or clum p should be primarily to distort it, rather 
than support it (|Ballesteros-Paredes et al.lll999al ). In highly 
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supersonic turbulence. I Vazquez-Semadeni et al] (|2008l ) have 
recently noted that the non-thermal motions within clumps 
must contain a significant compressive component, even if 
the driving is purely solenoidal. This compresive component, 
rather than opposing gravity, aids it or is driven by it. This 
is consistent with recent suggestions that, for example, the 
Orion A cloud is collapsing, an d producing the Orion Nebul a 
Cluster (ONC) in the process (|Hartmann fc Burkcrt 20071). 
that the non-thermal motions within the clump NGC 2264- 
C correspond mainly to gra vitational contraction, rather 
than to isotropic turbulence l|Peretto. Hennebelle fc Andrei 
l2007h . that massive star-formi ng regions may be immersed 
in large-scale accret i on fiows ([Ga lvan-Ma drid et al.l l2009l : 
ICsengeri erallboiol : [Schneider et al...2010. ), and that MCs 
in the Large Magellanic Cloud seem to follow an evolution- 
ary t rend such that m ore evolved clouds are more mas- 
sive (iFukui et al.ll20()9h . Also, if the non-thermal motions 
within clouds were homogeneous and isotropic turbulence, 
it would be difficult to understand the common observation 
that SF occurs at localized spots within the clo uds, rather 
than scattered throughout their volumes (e.g., iKirk et al] 
I2OO6; Evans et al. 2009). 

These findings all suggest that an important, and per- 
haps even dominant, component of the nonthermal motions 
observed in MCs and their substructure is actually converg- 
ing flows, which may be driven by gravity or by external 
compressions. In fact, a model in which the nonthermal 
motions in MCs are a gravit ationally-driven m ass cascade 
has been recently proposed bv lField et al.l l|2008l) . Moreover, 
the collision of converging flows has been shown to pro- 
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turbulent kinetic energy observed in objects as diverse as 
galactic disks, MCs, and protostellar accretion disks is in 
general consistent with being driven by infall from the en- 
vironments of those objects. However, if the turbulence is 
being driven by the gravitational contraction, it cannot be 
expected to halt the contraction that drives it. Thus, one is 
faced with the Zuckerman-Palmer conundrum again. 

One way to avoid the Zuckerman-Palmer conundrum is 
if large chunks of the molecular gas in the Galaxy are indeed 
magnetically subcritical and thus supported against gravity, 
while SF occurs preci sely in those regio ns that are not, as 
recently suggested bv lElmegreenI (|2007l ). Thus, the chaotic 
spatial and statistical distributions of the physical variables, 
produced by the turbulence in the forming MCs, may play 
a key role in the control of the SFR and the SFE. 

In this paper we present numerical simulations of the 
formation and evolution of MCs, starting from their for- 
mation out of generic compressions in the warm neutral 
medium (WNM), and reaching up to their star-forming 
epochs, in the presence of magnetic fields and AD, in or- 
der to investigate the production of sub- and supercritical 
regions, and the rate at which clouds with various environ- 
mental conditions form stars. We focus on the effect of AD 
on the SFE and the global evolution of the clouds. The plan 
of the paper is as follows; in sec. [2] we present some general 
considerations on the evolution of the M2FR in MCs. In sec. 
Owe present the numerical model and the parameters of the 
simulations. In sec. |4]we then present our results on the vari- 
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ability of the M2FR, the evolution of the SFE, and the role 
of AD in the evolution of both sub- and supercritical clouds. 
Finally, sec. [5] presents a summary and our conclusions. 



2 EVOLUTION OF THE M2FR IN 

MOLECULAR CLOUDS: FROM SUB- TO 
SUPERCRITICAL 

Although the most commonly considered mechanism for 
increasing the M2FR of a certain density enhancement is 
the redistribution of magnetic flux among the central flux 
tubes of a cloud by AD (jMestel fc Spitzeij|l956l : lMouschoviasl 
I1977I ). another important, yet often- neglected mechanism is 
that, for a uniform medium permeated by a given mean 
magnetic fleld strength Bq, there is always a certain length 
along the field (termed the "accumulation length" by Mestel 
1985; see also Shu et al. 2007) such that fiux tubes longer 
than that contain enough mass per unit area to be magnet- 
ically supercritical. The criticality condition in terms of the 
mass column density E = pL and the field strength Bp for 
a cylindrical geometry is jNakano fc Nakamuralll978l '). 



E/Bo « (47r^G) 



-1/2 



0.159 G" 



-1/2 



(1) 



where p is the mass density and L is the cylinder length. This 
condition gives the accumulation length, in terms of fiducial 
values representative of the ISM in the solar neighborhood, 
as |Hartmann ct al. 2001) 



' (5^0) Vl cm-3 



pc. 



(2) 



where n — p/{pmii) is the number density of the medium, 
mn is the Hydrogen mass, and is mean particle weight, 
taken as p — 1.27. In principle, if the Galactic field is pri- 
marily azimuthal, then the Galactic ISM at large is mag- 
netically supercritical in general, because field lines do not 
end, and thus sufficiently long distances are always avail- 
able along theml3 Thus, the M2FR of a system is not a 
uniquely defined, absolute parameter, but rather depends on 
where the system 's boundaries are drawn. We also stress that 
the M2FR depends on the local geometry of the consid- 
ered system. For instance, a system with spherical symme- 
try has a slightly lower critic al value of pcrit ~ 0.13 G"^'^'^ 
l|Mouschovias fc Spitzerlll976l ). Measuring the criticality of 
the streams with respect to this value would lead to a super- 
critical configuration for our runs B3 rather than sub-critical 
configurations. 

Now consider a cloud or clump that is formed by the 
accumulation of gas along field lines in general|f] In the rest 
of this discussion, we will generically refer to the result- 
ing density enhancement as a "cloud", referring to either 



^ Note, however, that supercriticality does not necessarily imply 
collapse, since the gas may be thermally or otherwise supported, 
as is likely the case for the diffuse warm medium at scales of 
hundreds of parsecs. 

^ Since compressions perpendicular to the magnetic field cannot 
induce collapse, and compressions oblique to the field can produce 
colla pse by reorienting the directions of the flow and the fleld lines 
llHen ncbcllc & Pcrault 2 000, ). our assumed configuration involves 
no loss of generality. 



a cloud, a clump, or a core. Although redistribution of mat- 
ter along field lines does not in principle affect the total 
M2FR along the full "length" of a fiux tube, this length 
is a rather meaningless notion, since the fiux tube may ex- 
tend out to arbitrarily long distances. What is more mean- 
ingful is the M2FR of the dense gas that makes up the 
cloud, since the cloud is denser than its surroundings, and 
thus it is the main source of the self-gravity that the field 
has to oppose. In fact, for the formation of a cloud out of 
flow collisions in the WNM, the clo ud's density is ~ 100 
times l arger than that of the WNM (jHennebelle & Peraum 

03; iHcitsch et al. 20(a 
VazQu cz-Semadcni et al.l 120061 : 
Baneriee et al. ,2009 ). and so the lat- 
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ter's self-gravity is negligible. Thus, in this problem, natural 
boundaries for the system are provided by the bounding sur- 
face of the dense gas, allowing a clear working deflnition of 
the M2FR. 

However, contrary to the very common assump- 
tion of a constant cloud mass, the formation of 
clouds by converging gas streams implies that the 
cloud's mass is a (generally inc r easing) function of time 
jBallesteros-Paredes et al. 1999bl: Vazguez-Semadeni et al.l 
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ceived observational support (jFukui et a l."2009'l. This means 
that, within the volume of the cloud, the M2FR is also an 
increasing quantity, since the flux remains constant if the 
flow is along fie ld lines, while the mass increases (see also 
IShu et al] |2007| ). If the cloud starts from essentially zero 
mass, this in turn implies that the M2FR of a cloud is ex- 
pected to start out strongly subcritical (when the cloud is 
only beginning to appear), and to evolve towards larger val- 
ues at later times. Rewriting eq. for the column density, 
we see that the cloud becomes supercritical when 



(3) 



where A'^ = "E/pmH is the number column density along 
field lines. 

The critical column density for magnetic criticality 
given by eq. @ turns out to be very similar, at least for 
solar neighbourhood conditions, to the critical column den- 
sity of hydrogen atoms necessary for cold at omic gas to be- 
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that it starts out as an atomic and subrcritical diffuse cloud 
i Vazquez- Semadeni et ai\\200(i) and, as it continues to ac- 
crete mass from the warm atomic medium, it later be- 
comes molecular and s upercritical, roughly at the same time 
(|Hartmann et al]|200ll ). This is fully consistent with the ob- 
servation t hat diffuse atomic cloud s are in general strongly 
subcritical (|Heiles fc Trolandll2005l '). while MCs a rc approx- 
imately critical or moderately supercritical (jCrut chcr 199^; 
iBourke et allbOOll ; FlVoiand fc Crutcheill2008D . 

Moreover, the critical column density given by eq. ([HJ 
is also very similar to that required for rendering cold gas 
gravitationally unstable, which is estimated to be 



4 Vazquez- Semadeni et al. 



0.7 X 10^ 
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3000 K cm-3 



1/2 



(4) 



l|Franco fc Coxlll986l : [Hartmann eral] 120011 ). Thus, at so- 
lar neighborhood conditions, a forming cloud is expected 
to become molecular, magnetically supercritical, and self- 
gravitating at roughly the same time. 

It is important to note that the mass accretion onto 
a cloud due to gas stream collisions is likely to start 
along essentially just one dimension. This mode of mass 
accretion may be slow, and it has be en argued that it 
may involve excessively long times (e.g.. iMcKee fc Ostrikeii 
I2OO7I ). However, numerical simulations of the process show 
that, once the gas has transitioned to the cold, dense 
phase, it soon becomes gravitationally unstable, even 
though it may remain mainly in the atomic phase, and 
three-dimensional gravitational contraction can then en- 
sue, providing a much faster mode for inc r easing the 
column densit y ('Vazguez- Semadeni et af] 20071: ElmegreenI 
20071: iHeitsch fc HartmannI I2OO8I : iHennebelle etal.1 I2OO8I : 
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vergence of the flow is driven by large r -scale gravita- 
tional instabilit i es (e .g., iKim et all l2003l : iLi et all l2005l : 
iKim fc Ostrite3[2007l '). Of course, this increase of the col- 
umn density due to gravitational contraction of the dense 
gas is only relevant to molecule formation. During such a 
process, the M2FR remains constant if the cloud's mass re- 
mains fixed or varies on timescales much longer than the 
contraction, and the latter occurs under ideal MHD condi- 
tions. The gravitational contraction can only contribute to a 
further increase of the M2FR if the gravitational potential of 
the cloud causes it to accrete further amounts of diffuse gas, 
which transitions to the dense phase as it is incorporated 
into the bulk of the cloud. 

In all of the processes discussed so far, AD has not 
played a role. This is of course due to the well known fact 
that AD is not dynamically relev ant until densities as hig h 
as nAD ~ 10^ cm~^ are reached (|Mouschovias et al.lll985r ). 
Such densities are only reached in the dense cores of MCs 
and therefore AD is not expected to cause any important 
effects on the global evolution of MCs. Such cores, however, 
may be magnetically subcritical if they form by turbulent 
compressions within the cloud before AD becomes locally 
important, even if the cloud is globally supercritical. This is 
because, under ideal MHD, a core formed within an initially 
unifor m cloud must have a smaller M 2FR than that of the 
cloud (|Vazauez-Semadeni et af]|2005l ). In a sense, the core 
repeats the pattern followed by its parent cloud, initially be- 
ing strongly subcritical and evolving towards higher values 
of the M2FR, being limited by the M2FR of its parent cloud, 
until AD becomes important and allows its M2FR to over- 
take that of the parent cloud, perhaps becoming supercrit- 
ical and allowing the core to collapse. However, this notion 
has not been tested in the context of the global evolution of 
a GMC, in particular taking into account the property that 
the cloud's M2FR should evolve (generally increasing) with 
time. In the remainder of the paper we investigate this sce- 
nario, by means of numerical simulations of the formation 
and evolution of a GMC, including AD, and focusing in par- 
ticular on the resulting SFE. We are particularly interested 
in the star-forming properties of the cloud as it transits from 
sub- to supercritical. 



3 THE NUMERICAL MODEL 

3.1 The numerical code and setup 

We use the adaptive mesh refinement (AMR) code FLASH 
(|Frvxell et al.l bood ) wit h MHD, modified to in clude the 
AD module developed bv lDufhn fc Pudritj l|2008l ). The AD 
treatment takes the single-fiuid approximation, and uses a 
simple prescription to avoid the need to track the ion density 
in this approximation. This prescription essentially turns off 
AD at low (n < 1 0'^ cm~^) densitie s . For more details, we 
refer the reader to lDuffin fc Pudritj l|2008l ). 

The simulat i ons also use a sink particle prescription 
llBate et al.lll995l; iJappsen et al.|[2005l : iBaneriee et"al]|2009l : 
iFederrath et al.l 2010l ). A sink particle is created in a cell 
if the density there reaches a threshold density risink ~ 
2 X 10^ cm~^, and the cell is a local minimum of the gravi- 
tational potential. When a cell forms a sink, the latter takes 
all the mass in excess of Tisink in the region where the density 
n satisfies n ^ risink- The sink particles have an accretion 
radius of 0.065 pc, corresponding to roughly 1 Jeans length 
at risink and T ~ 20 K. 

Concerning the heating and cooling, we u s e the same 
prescription us e d in IVazguez-Semadeni et al.l (|2007l ) and 



I3aneriee et al. [20(M) , wh ich is derived from the fit by 
Kovama fc lrmtsukaT (l2002|') to the results of the che mistry 
and cooling calculations of Kovama fc Inutsukal l|2000l l . This 
prescription implies that the simulated ISM is thermally un- 
stable in the density range 1 J; ^ 10 cm^"^, which, under 
thermal balance between heating and cooling, corresponds 
to the temperature range 5000 > T > 500 K. 

We model the convergence of WNM flows as the colli- 
sion of two large-scale cylindrical streams. Our setup is sim- 
ilar (t hough not identical) to the non-m agnetic SPH simula- 
tion of IVazquez-Semadeni et all |20o3) labeled L256Ar;0.17 
(see their Fig. 1). Each stream is 112 pc long and has a ra- 
dius of 32 pc. The streams collide at the plane x = pc, 
and are embedded in a (256 pc)^ simulation box, in which 
the coordinates range from —128 to 128 pc. The numerical 
box is periodic, and the streams are completely contained 
within it, ending at a distance of 16 pc from the x bound- 
aries. The resulting cloud occupies a relatively small volume 
far from the boundaries, and so it can interact freely with 
its diffuse environment, with relatively little effect from the 
boundaries. Most importantly, the cloud is free to grow by 
accretion from the WNM. 

The cylindrical streams are given an initial, moder- 
ately supersonic inflow velocity Vinf so that they collide at 
the centre of the numerical box. The inflow speed of each 
stream is measured with respect to the isothermal sound 
speed of 5.7 km s~^ that corresponds to the initial temper- 
ature of 5000 K. This isothermal inflow Mach number is 
denoted A^inf. We also add 10% random velocity pertur- 
bations to the bulk stream speeds, in order to trigger the 
instabilities tha t generate turbulence i n the forming cloud 
(IVishniad Il994l: iHeitsch et al.l l2005l : iPittard et ail l2005l : 
IVazquez-Semadeni et alT 20061 ). The box is initially filled 
with WNM at a uniform density of n = 1 cm~^ (p — 2.12 x 
lO^^** g cm~^, using a mean atomic weight of 1.27). At the 
temperature of T = 5000 K for the warm phase, this implies 
that the cold phase comes into hydrostatic thermal pressure 
bala nce with the WNM at a density n ~ 100 cm~^ (see Fig. 
2 of IVazguez-Semadeni et al]|2007h . However, in our simu- 
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lations the density of the cold phase is higher, because it is 
in balance with the sum of the thermal and the ram pres- 
sure of the col liding streams ( Vazgucz-Sc madeni et al ] |2006l : 
lHennebelle"etal.,.2008, : ,Baneriee et al...200^ 



sion will again be to allow SF to occur more readily than if 
mediated by AD alone in subcritical cases, and so our SFRs 
must again be considered upper limits to the ones caused by 
AD alone. 



3.2 Resolution issues 

We start our simulations at a base resolution of 512^, i.e. 
Ax = 0.5 pc, at the convergence point of the flows. Addi- 
tonally, we allow the code to refine up to 4 additional levels, 
the highest of which corresponds to a maximum resolution 
of 8192 grid points, or a grid spacing of Ax = 0.03 pc in 
each direction. For t he dynamical mesh r efinement we use a 
Jeans-type criterion l|Truelove et al ] |l997l . see, however, Fed- 
errath et al. (2011, submitted) for a more stringent criterion 
in the presence of magnetic fields) , requiring the local Jeans 
length to be resolved with at least 10 grid cells while refining 
is active. Beyond the last refinement level, the Jeans length 
begins to be more poorly resolved, until a maximum allowed 
density of risink = 2 x 10^ cm~^ is reached, at which a sink 
partic le is formed. We refer the reader to iBaneriee et al.l 
(|2009l ) for a discussion of the justification and possible lim- 
itations of this choice with regards to thermal issues. Here 
we discuss issues related with gravity and AD. 

The value of risink we use was chosen in order to reason- 
ably ensure that the material going in to sink particles is ac- 
tually gravitationally bound. Indeed, ICalvan-Madrid et al.l 
(|2007l ) found that, when cores are defined by means of 
a density threshold, most cores defined by thresholds > 
10^ cm~^ proceeed to collapse. Instead, when cores are de- 
fined by lower thresholds (say, ~ 10* cm~^), a significant 
fraction of th em is transient, rebounding instead of collaps- 



ing ( see also iBallesteros-Paredes et al ] |2003l :lKi essen et al.l 
I2OO5I : IVazauez-Semadeni et all I2OO5I ) 



Now, at n — risink 
and T = 20 K, the Jeans length is Lj = 0.066 pc, and 
so we marginally fail to fulfill the minimum Jeans criterion, 
of resolving Lj with at least 4 cells. However, this should not 
introduce any significant errors, as we are not concerned here 
with the fragmentation of the core into multiple stars, nor 
with their mass distribution, but only with the total mass 
going into stars. 

Another issue is that numerical diffusio n can have an 
effect similar to that of AD, as discussed bv iKlessen et al.] 
l|2000h . Specifically, since the scale of the densest cores is a 
few grid cells, numerical diffusion can cause the magnetic 
flux to diffuse out of them, in a similar manner to AD. In- 
deed, we do occasionally observe the occurrence of gravita- 
tional collapse in magnetically subcritical, ideal MHD simu- 
lations, in which theoretically this should not occur. More- 
over, since risink effectively constitutes an upper limit to the 
density that can be reached by any cell in the simulation, 
and since we have chosen a value of risink that is of the same 
order as riAD, the highest densities in the code will be of 
the order of riAo, and numerical diffusion and AD will have 
comparable effects. This limitation could be avoided by us- 
ing an even larger number of refinement levels but, since 
the simulations are already very numerically expensive (~ 
200,000 CPU hours per run), this option is not presently 
feasible. Alternatively, we could give up on satisfying the 
Jeans criterion, simply raising risink without increasing the 
allowed number of refinement levels, an option that we may 
attempt elsewhere. In any case, the effect of numerical diffu- 



3.3 The simulations 

We consider five numerical simulations with three reason- 
ably realistic values of the initial, uniform magnetic field 
-Bo, of 2, 3, and 4 /iG, respectively. These values span the 
observed range of valu es of the un iform component of the 
Galactic mag netic field (|Beckll200ll ). The initial field is along 
the x-direction. 

With respect to the cylindrical criticality criterion, eq. 
lU, these cases respectively correspond to fi/ficrit ~ 1.36, 
0.91 and 0.68, so that the first case is magnetically super- 
critical while the other two are subcritical. Note, however, 
that the subcritical cases are only so because of the finite 
extent (256 pc) of the numerical box. For Bo = 3 and 4^G, 
lengths of 280 and 380 pc, respectively, would be required 
to render the system magnetically critical. Also, because the 
critical value of the mass-to-flux ratio depends on the geom- 
etry of the considered configuration, individual (molecular) 
clumps could be supercritical if compared to the slightly 
lower ^crit of iMouschovias fc Spitzeii (|l976l ) (see discussion 
in SecHi}. 

The supercritical case is considered only in the AD 
regime, as we do not expect the absence of AD to make 
a significant difference in this case. The _Bo = 3 and 4/xG 
cases are considered both in the "ideal" and AD regimes, to 
investigate the effect of AD on the star-forming properties of 
magnetically subcritical clouds. Note that we have written 
the word "ideal" within quotation marks because we cannot 
avoid the effect of numerical diffusion, even if we turn off 
the AD. Except for the value of the magnetic field strength 
and whether AD is on or off, the simulations are otherwise 
identical, all having an inflow speed of 13.9 km s~^, corre- 
sponding to an isothermal Mach number A^inf = 2.44 with 
respect to the unperturbed, initial medium at T = 5000 K. 

The runs are labeled mnemonically, so that the first two 
characters of the run's name indicate the field strength in 
piG (e.g., "B3" denotes Bq = 3/iG), and the last two denote 
whether the simulation is in the ideal MHD case ("MH") 
or includes AD ("AD"). Table [1] summarizes the parameters 
used in each of the five runs. The last column in this Table 
gives the maximum time reached by each simulation. 



3.4 Considerations on measuring the mass-to-flux 
ratio 

In what follows, we will be presenting measurements of the 
M2FR in various regions of the simulations. However, this is 
not an unambiguous task in general, and in fact the M2FR 
can be measured using different procedures. In principle, the 
M2FR should be measured along flux tubes, in order for the 
measurement to be directly representative of the dynamical 
effect of the field on the gas. Thus, the measurement should 
be performed by tagging a bundle of field lines, and inte- 
grating the density along the path defined by them. Unfor- 
tunately, such a measurement is extremely difficult to per- 
form, even in the simulations. A magnetic flux tube may lose 
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Table 1. Run parameters 



Run Bo AD M/A'crit Final time 
name [fiG] [Myr] 



B2-AD 


2.0 


on 


1.36 


31.4 


B3-MH 


3.0 


off 


0.91 


26.1 


B3-AD 


3.0 


on 


0.91 


35.6 


B4-MH 


4.0 


off 


0.68 


48.4 


B4-AD 


4.0 


on 


0.68 


59.2 



Table 2. Initial conditions and the final simulatation time of the 
runs presented in this work. To calculate the magnetic critically 
of the entire system we use the critical mas s-to-flux ratio, /^crit ~ 
0.16G ~^/^ for a cylindrical geometry of iNakano Sz Nakamural 
b97^ 



coherence if the field lines that compose it diverge from each 
other at long distances. Also, the field near and within the 
cloud can be significantly distorted, due both to the turbu- 
lence in the cloud, and to its gravitational contraction, even 
if the initial flow direction is along the field lines. 

Observationally, the M2FR is often estimated by mea- 
suring the ratio of column density to magnetic field strength, 
N/B, alon g lines of sight (LOSs ) through the cloud of in- 
terest (e.g. ICrutcher et al.ll20o3 ). However, this procedure 
actually intersects many different fiux tubes, and thus gives 
only an approximation to the ac t ual M 2FR of a single flux 
tube. As discussed by ICrutcheil () 19991 ). if a cloud is flat- 
tened, its plane is perpendicular to the magnetic fleld lines, 
and the system is observed at an angle 0, then the observed 
M2FR, or equivalently, the N/B ratio, will be related to the 
actual one by {N/B)ohB ~ N/{B cos^ 6). Alternatively, as is 
the case for the measurements we present below, if the ob- 
servation is performed along an LOS that is perpendicular 
to the cloud, but the magnetic fleld is at an angle 9 with 
respect to the LOS and to the normal to the plane, then we 
have that (iV/B)obs = iV/(B cos 61), with l^l < 7r/2. In ei- 
ther case, {N/B)oi,s tends to overestimate the actual value, 
and on some occasions very large values may be artiflcially 
measured. This implies that a map of {N/B)obB i s actually 



a map of upper limits to the actual N/B. This led lCrutcheij 
(|l999l ) to introduce statistical correction factors of 1/2-1/3 
to the set of M2FR values obatined in the observations he 
considered. 

The equivalent procedure for the simulation data is to 
measure M2FR along LOSs through the clouds in our sim- 
ulations. We choose the LOSs to lie along the a;-direction, 
since this is the direction of the mean magnetic fleld and 
of the colliding WNM streams, and thus it is the direc- 
tion along which the column density is dynamically relevant. 
Specifically, we then measure the M2FR as 



(5) 



where L is the stretch of the cloud along the x direction. 
The path L is chosen so as to contain the full extent of 
the cloud's thickness. This is done for every position on the 
plane of the cloud, to obtain maps of the M2FR over the 
cloud's surface. We refer to this as "the projection method" 
of measuring the M2FR in the simulations. 



In order to appreciate the amount of distortion that may 
be present in the field within the clouds, in Fig. [1] we show 
cross sections through the centre of runs B3-MII and B4- 
MH along the (x, y) plane, showing the density field and the 
component of the magnetic field on this plane. We observe 
that in run B3-MH the field lines are not strongly bent, 
except at the sites of local collapse. This suggests that the 
M2FRs we measure by the projection method should not 
exceedingly overestimate the actual fiux-tube value. 

Nevertheless, a more definite way to estimate the 
amount of overestimation of the M2FR incurred in by the 
projection method is to use a different estimator. One such 
estimator is what we call the "local method" , which consists 
in measuring the M2FR for individual cells in the simula- 
tions, using the total magnetic field strength. Specifically, 
we measure the gas mass M and read off the total magnetic 
field strength _B in a cell, in order to calculate the M2FR as 

^^^Bd^' 

where dx is the cell's side length. This estimator is gives a 
lower bound to the M2FR in a magnetic fiux tube, since it 
only counts the mass in a single cell within that tube. This 
method has no observational analogue but, by using the two 
methods, we expect to bracket the true distribution of values 
of the M2FR, at least in a statistical sense. 



4 RESULTS 

4.1 Global evolution and star formation 

We first direct our attention to the global evolu- 
tion of the clouds. The supercritical run B2-AD 
evolves very s imilarly to the non-magnetic 
presented by 



Vazguez-Semadeni et al.l 

IVazguez-Semadeni et al.l 1 20101) and the 
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and 

su- 

percritical run s prese nted by Hennebelle etal.1 l|2008l ) and 
iBaneriee et al.l (|2009l ). The cloud starts out as a thin 
cylindrical sheet that fragments and thickens as time 
increases, until it becomes gravitationally unstable and 
begins a global radial contraction at t ~ 9 Myr. Shortly 
after that {t ~ 12 Myr) star formation begins in the densest 
fragments ("clumps"), while the fragments continue to fall 
towards the global centre of mass, and by t 24 Myr a 
dense cloud of radius ~ 10 pc has formed there, which does 
not appear to contract further. This lack of contraction, 
however, is only apparent, because in fact gas is being 
consumed within the cloud by SF, and gas from the outside 
continues to fall onto the cloud. Figure [2] shows this run at 
t = 10, 20 and 30 Myr, illustrating its evolution. 

On the other hand, the subcritical runs B3 and B4 un- 
dergo a period of initial contraction followed by a rebound, 
eventually settling into an oscillatory regime, which consists 
of alternating periods of contraction and expansion around 
th e magnetostatic equilib rium state, as previously observed 
bv lLi fc Nakamural (|2004l ) . These oscillations are best seen 
in animations of the simulations (not shown), but they can 
also be observed in Figs. [3]and|4l which show snapshots of 
the density field of runs AD-B3 and AD-B4 at various times, 
respectively. In Fig. [31 it can be seen that the central density 
is larger at the intermediate time shown in the middle panel 
than at the flnal time shown at the right panel. A similar 
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behavior is seen in Fig. |4l where the central density is seen 
to be larger at times t = 20.5 and i = 48 Myr than at t = 34 
Myr. The B3 runs, having a weaker mean field, contract for 
a longer time (up to f ---^ 25 Myr) and reach a smaller size 
(radius i? ~ 20 pc) than the B4 runs (maximum contraction 
at t ~ 20 Myr, with radius 7? ~ 25 pc). The B4 runs, which 
were followed to longer times, clearly exhibit the oscillatory 
regime, with a period of ~ 30 Myr. The oscillations can also 
be seen in Fig. [S] which we now discuss. 

Figure [S] shows, for all the runs, the evolution of the 
total dense gas mass and total sink mass (top left panel), 
the time derivative of the total sink mass Msinks (top right 
panel), the total number of sinks {bottom left panel), and the 
SFE {bottom right panel), defined as 



SFE(t) = 



Meinks(t) 



A7dcnsc(i) + Msinks{t) ' 



(7) 



where Mdcnsc is the mass of the gas with density n > 
100 cm~^, and Mainks is the total mass in sink particles. 
We take Mginks as a proxy for the SFR. 

The evolution of the runs is seen to depend sensitively 
on the magnetic field strength. For t > 7 Myr, the dense 
gas mass oscillates by factors of 2-4 for the subcritical runs 
B3 and B4, in both the MHD and the AD cases, evidencing 
again the oscillatory regime in which these runs engage. The 
maxima of the mass in the subcritical runs coincide with the 
times of maximum compression. This may be partially an 
artifact of the threshold density we have chosen for defin- 
ing the cold gas. Because the global magnetic support for 
the cloud prevents it from contracting much, the gravita- 
tional binding of the cloud is generally weak. This in turn 
means that the during periods of maximum expansion, a 
significant fraction of the cloud, although still in the cold 
phase, may be below the density threshold of 100 cm~^ we 
have used for defining the cloud. Nevertheless, if the mean 
density of the cloud varies, it is likely that the molecular 
fraction should actually vary as well, since molecular gas 
may be diss ociated if the clou d's column density decreases 
sufficiently (iGlover et aLlbOld ). Thus, during periods of ex- 
pansion, the cloud may contain a lower molecular fraction. 
In addition, significant amounts of gas may be in a transient 
state betw een the w arm an d cold phases of the ISM (see the 
reviews bv I Vazque z- Scma deni et al.ll2003l : iHennebelle et all 
I2OO9I : IVazauez-Semadenill20'09r and references therein), be- 
ing thus lost from the cold, dense phase. 

The supercritical run B2, on the other hand, does not 
exhibit such strong oscillations in its dense gas mass content. 
Instead, the dense gas mass increases rapidly at first. This is 
because this run does not engage in any radial oscillations, 
but simply proceeds directly to collapse. Interestingly, how- 
ever, the dense gas mass later becomes roughly stationary, 
although this stabilization is not due to the cloud being in 
any sort of equilibrium, but rather to the fact that it is form- 
ing stars at roughly the same rate it accretes mass from the 
WNM. Indeed, the red lines in the top left panel, as well as 
the top right panel of Fig. [S] show that the total sink mass in 
the simulation increases at a steady pace, of roughly 400Mq 
Myr"^ during the time interval 20 < t < 30 Myr. 

All runs, including those that are subcritical, form 
"stars" (i.e., sink particles), since, as discussed in Sec. 13.21 
numerical diffusion acts in a similar manner to AD. Never- 
theless, when AD is included the total sink number and mass 



increase, the dense gas mass decreases, these effects being 
relatively more noticeable for the strong-field runs B4. In 
general, as can be seen from the bottom left panel of Fig. [S] 
the total number of sinks in the B4 and B3 cases is larger 
only by a few sinks when AD is included. This reinforces 
our conclusion from Sec. l3.2l that numerical diffusion has an 
effect of comparable strength to that of AD in our simula- 
tions, since AD is able to induce the collapse of a few clumps 
in addition to those that collapse because of numerical dif- 
fusion. In the case of the B3 runs, the relative effect of AD 
is smaller because the same difference of a few extra sink 
particles is a smaller fraction of the total number of sinks 
produced. 

Indeed, run B3-AD forms stars at a much higher rate 
than run B4-AD, even though both are subcritical. The time 
derivative of the sink mass, i\fsinks, in run B3-AD is roughly 
one order of magnitude larger than that of run B4-AD, as 
seen in the top right panel of Fig. [S] In fact, the formation of 
sinks completely stops in run B4-AD after t ~ 27, as can be 
seen in the bottom left panel of this figure. The very slight 
increase in the sink mass observed during this period {top 
left panel, dash-dotted red line) is due to accretion onto the 
existing sink particles, rather than to the formation of new 
particles. Run B2-AD, on the other hand, is seen to have a 
larger SFR than run B3-AD, although only by factors of a 
few. Also, it can be seen from all panels of Fig. [5] that the 
onset of SF is delayed as Bo is increases, but that the pres- 
ence of AD shortens this delay. It is important to note that 
all SF in the runs occurs after the cloud has been assembled 
and its M2FR has reached a nearly stationary value {t > 10 
Myr). 

To conclude this section, the bottom right panel of Fig. [5] 
shows the evolution of the SFE, defined by eq. , in all sim- 
ulations. Again, a continuous trend of increasing SFE with 
decreasing Bo is observed throughout our set of simulations. 
The supercritical run B2-AD reaches an SFE of 35% at 
t = 30 Myr. For comparison, at this time, run B3-AD has 
reached an SFE of ~ 25%, while run B4-AD has reached only 
SFE ~ 3%. Although at face value these numbers would sug- 
gest that run B4-AD compares best to the ob served SFEs 
of GMCs (|Mvers et al.ll 19861 : lEvans et al.ll2009l '). this conclu- 
sion may be premature, since the additional effect of stellar 
feedb ack in reducing the SFE (e.g.. lVazauez-Semadeni et al.l 
is not taken into account in the present simulations. 
We conclude that the SFR and the SFE can depend sensi- 
tively on the mean field strength, even for globally subcrit- 
ical cases, and that the marginally subcritical run has an 
SFE comparable to that of the supercritical case. 



4.2 Spatial and probability distribution of the 
mass-to-flux ratio 

A key piece of information needed to understand the behav- 
ior of simulations is the spatial distribution of the M2FR, 
as well as the evolution of its global average value. In what 
follows, we discuss the M2FR estimated using the projection 
method (cf. eq. [S]) along LOSs parallel to the x axis (i.e., 
perpendicular to the plane of the cloud), taking L as the 
path — 10 < x < 10 pc along the direction of the infiows. We 
consider the M2FR normalized to the critical value given 
by eq. ([T]). Given the flattened geometry of our clouds, we 
consider that this is a more realistic value of the critical 
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M2FR than the other commonly encountered critical value 
of (67r^G)"^/^ « al3G"^/^ which holds for spherical ge- 
ometry (e.g., 'Shul ll992l ). Figure [6] shows snapshots of the 
normalized M2FR for runs B2-AD, B3-AD and B4-AD in 
the top row, and of runs B3-MH[ and B4-MH in the bottom 
row, all at t = 20 Myr. In all cases, the spatial distribution 
of the M2FR is seen to fluctuate strongly. 

Comparing the MH and AD cases, it is interesting to 
note that the spatial structure of the M2FR is similar in 
the large-scale features, but differs in the shape and precise 
location of the fine, small-scale ones. This is a manifestation 
of the system being chaotic, so that the subtle variations 
in the magnetic forces at the densest structures induced by 
AD are sufficient to change the details in the topology of the 
gas. Presumably, for sufficiently long times the differences 
in structure will reach even the largest scales. On the other 
hand, simple visual inspection of the images is not enough 
to discern any trend of systematically larger values of the 
M2FR in the presence of AD. To quantify this, we show 
in Fig. [7] the histograms of the M2FR. The histograms are 
computed for all LOSs within a circular region centered at 
{y,z) — (0,0), with a 20-pc radius. From these, we see that 
the inclusion of AD causes the production of a small excess 
of high-M2FR cells in comparison with the non-AD cases, 
and that this effect is most noticeable in the strong-field 
(B4) case, in which the maximum value of the M2FR is over 
a factor of 2 larger than in the non-AD case. Instead, in 
the B3 case, the excess is marginal, suggesting that when 
the system is very close to being supercritical, numerical 
diffusion dominates over AD. This result suggests that the 
relative importance of AD and numerical diffusion depends 
on the mean field strength, an issue that deserves further 
exploration, but which is out of the scope of the present 
paper. 

Returning to Fig. [G] and focusing on the top row of 
images, which show the variation of the M2FR's spatial 
distribution as a function of the magnetic field strength, 
several points are worth noting. First, as mentioned above, 
the M2FR is seen to be highly inhomogeneous in all three 
runs. This is also illustrated in Fig. [S] whose top left panel 
shows histograms of the M2FR in the three runs at the 
same time as that shown in Fig. (5] The top right and bot- 
tom left panels of Fig. [S] show the corresponding cumu- 
lative distributions, respectively weighted by volume and 
by mass. Finally, the bottom right panel shows the mass- 
weighted cumulative distribution for high-column density 
{N > 10^^ cm"^) LOSs only. From these fi gures, it is seen 
that fj, fiuctuates by over one order of magnitude in the 
subcritical runs, and by two in the supercritical one. Part 
of this variability, especially the highest values of fi, may 
be an artifact of the measurement procedure, as discussed 
in Sec. 13.41 Nevertheless, significant actual fluctuations of 
the M2FR on the plane are expected, since the cloud is 
turbulent and clumpy, due to the combined action of the 
thermal, Kelv in-Helmholz a nd nonlinear thin-shell (NTSI , 
Vishniag 19941) inst abilities (|Vazauez-Semadeni et af]|2006l : 
Heitsch et al.1 120061) . In the ideal MHD case, segments of 
magnetic flux tubes must have M2FRs smaller than that 
of the whole tube at all times. However, in the presence of 
diffusion (numerical and/or ambipolar), Lagrangian regions 
where the density is large enough can lose magnetic flux 
and reach M2FR values larger than the initial value for the 



tube. Thus, local clumps may become magnetically super- 
critical even within the globally subcritical simulations, as 
in the low-mass m ode of the "standar d model" o f magneti- 
cally regulated SF (|Shu et a l.|| 1987: Mouschoviasi,199 1). and 
turbulent extensions of it i Nakamura fc Lil l2005l )T Accord- 
ing to the top right and and bottom left panels of Fig. [S] the 
fraction of the volume (resp. mass) that is in locally super- 
critical regions increases from ^ 3% (resp. ~ 10%) in run 
B4-AD to ~ 17% (resp. - 58%) in run B2-AD. In summary, 
there exist plenty of mechanisms that contribute to the de- 
velopment of a highly inhomogeneous spatial distribution of 
the M2FR, both in the ideal MHD and in the diffusive cases. 
This is consistent with recent observational determinations 
suggesting that the magnetic field strength is randomly dis- 
tributed in MCs, with onl y its maximum value s scaling as a 
power law of the density (ICrutcher et al.|[2010l ). 

Second, we note from the top left panel of Fig. [S]that the 
M2FR histogram for the supercritical run B2-AD is wider 
than the histograms of either of the subcritical runs, having 
a larger fraction of both sub- and supercritical LOSs. This 
can be understood as a consequence of the combined ac- 
tion of AD enhanced by turbulence and mass conservation. 
For weaker magnetic fields, the Alfvenic Mach number is 
larger, implyi ng larger density fiuctu a tions, in which AD 
is enhanced (iFatuzzo fc Adamsl |2002| : iHeitsch et al] |2004| : 
iLi fc Nakamura ''2004')." thus allowing the formation of more 
strongly supercritical clumps, which are denser and contain 
more mass. In turn, this implies a stronger evacuation of the 
remaining regions, which are those left with smaller masses, 
and therefore with lower values of the M2FR. 

Third, from the images in Fig.[51 in which the dots indi- 
cate the positions of the sink particles, we see that not all of 
the regions that appear supercritical according to the pro- 
jection method proceed to form stars. This may be either be- 
cause they are truly supercritical albeit locally Jeans-stable, 
or because they are actually magnetically subcritical, and 
only appear supercritical due to the projection effect. There- 
fore, it is important to determine the degree to which the 
M2FR is overestimated by the projection method. We defer 
a detailed energy-balance study of the clumps and cores for 
a future paper, but here we can take a first step towards 
addressing this problem by comparing the maps and his- 
tograms of the M2FR obtained with the projection method 
to those obtained with the "local method" (cf. Sec. l3.4|) . Fig- 
ure |9] shows histograms of the M2FR using this method for 
the middle plane of each simulation at t = 20 Myr. Here we 
do not integrate over any LOS in order to show the largest 
excursions that the M2FR can exhibit in local cells. As men- 
tioned in Sec. 13.41 the local method gives lower limits to the 
actual M2FR in the fiux tube to which the local cell belongs. 

The /x-histograms using the local method exhibit a imm- 
ber of interesting features. First, it is seen that no super- 
critical cells are seen in both of the subcritical runs. At 
least for run B3-AD this is necessarily an underestimation 
of the actual M2FR, since sink formation has already oc- 
curred in this run at the time at which the histograms are 
made {t = 20 Myr). Second, the supercritical run B2-AD ex- 
hibits a small but finite fraction of supercritical cells, even 
with this M2FR-underestimating method, indicative of the 
abundance of supercritical regions in this case. Third, the 
histograms are seen to peak at /x/^crit ^ , while those 
obtained with the projection method peak at /i//icrit ~ 0.6. 
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Thus, the difference between the two methods is so large 
that the true distribution remains relatively unconstrained 
between them. 

A final criterion that can be used to determine the 
accuracy of the observational-like projection method is to 
compare the supercritical mass fraction obtained with this 
method with the SFE, which is in fact a measure of the mass 
fraction that has become simultaneously Jeans-unstable and 
magnetically supercritical over the cloud's history. Unfortu- 
nately, in principle this relationship is not trivial, because 
the supercritical mass fraction is an instantaneous quantity 
in the simulation, while the stellar mass in the cloud is a 
time-integrated quantity. However, if an approximately sta- 
tionary state is established in which fresh gas is continuously 
added to the clouds by accretion f rom the WNM, replenish- 
ing t he gas used up to form stars (jVazgue z-Semad eni et all 
l20ld ). then the (stationary) SFE may be considered as a 
lower limit to the (stationary) supercritical mass fraction. 
The notion that the SFE is a lower limit of the supercriti- 
cal mass fraction allows for the possibility that efhciency of 
conversion of gas to stars is still less than unity even within 
the supercritical, Jeans-unstable gas. 

With this in mind, we plot in Fig.[TO]the SFE versus the 
supercritical mass fraction in the three runs B2-AD, B3-AD 
and B4-AD, as read off from Figs. [5] {bottom right panel) 
and [8] {bottom left panel). The plotted value of the SFE is 
the mean between the extremes taken by the SFE over the 
time interval after which the initial rapid growth has ended, 
and the error bars denote these extremes. The dotted line 
indicates a least squares fit to the data points, given by 



SFE -0.034 + 0.54 (^^^^ . (8) 

Of course, this fit is totally empirical, and is provided only 
as a guideline for the trend of the SFE with the supercritical 
mass fraction. In particular, it is meaningless below the value 
of the supercritical mass fraction that produces an SFE of 
zero. 

We observe that, in all three cases, the SFE is a few to 
several times smaller than the supercritical mass fraction. 
Under the assumption of stationarity, this then suggests that 
the supercritical mass fraction obtained through the projec- 
tion method does not differ from the real one by more than 
factors of a few, on average. 

A final point worth noting is that, as shown in the bot- 
tom right panel of Fig. [SI the supercritical mass fraction 
in the high-column density gas is significantly larger in all 
three runs than the supercritical mass fraction for gas at all 
column densities {bottom left panel). This reinforces the sce- 
nario that the M2FR is determined mainly by the column 
density acquired by the individual regions in the clouds by 
accretion of gas along field lines, and less importantly by the 
local effect of ambipolar (and/or numerical) diffusion. 

We conclude from this section that the projection 
method gives estimates of the M2FR that are within less 
than a factor of a few from the actual distribution, and that 
the dominant mechanism that determines the local M2FR 
is the accumulation of gas along field lines. 



4.3 Evolution of the global M2FR 

We proceed now to discuss the evolution of the mean M2FR 
of the clouds in the simulations, and its range of variabil- 
ity. Figure [TT] shows the evolution of the mean and the 3a 
range of the normalized M2FR, /i//icrit, for runs B2-AD {top 
panels), B3-AD {middle panels), and B4-AD (bottom pan- 
els). The computation of the M2FR is performed using the 
projection method with the same path length and over the 
same circular region as those used for the histograms of Fig. 
[S] The left panels of Fig. [TT] show the mean and 3a range 
computed for the set of all lines of sight contained in the 
circular region, while the right panels show these quantities 
computed only for the set of lines of sight for which the 
column density is larger than lO'^^ cm~^. 

These plots illustrate the fact, discussed in Sec.[2l that 
the M2FR of a cloud is an evolving quantity, which first 
increases as the cloud gathers material from the WNM in- 
flows that assemble it. At the inflow speed of 13.9 km s~^, 
the 112-pc-long inflows are entirely incorporated into the 
cloud in approximately 8 Myr, which indeed corresponds to 
the time scale at which the M2FR of all three clouds is seen 
to have reached a roughly stationary value. This value is 
seen to be larger for weaker mean held, as expected. How- 
ever, a second increase in the M2FR is seen to occur after 
this first stabilization. This can be attributed to the fact 
that the inflows drag part of the surrounding, initially static 
WNM along with them, as they leave a partial vacuum be- 
hind them. This dragged material flows at lower speeds, and 
reaches the cloud at later tim es (for further discussion, see 
IVazguez-Semadeni et al.ll2007l '). 

It is interesting that {fi/fiait) for the set of all lines of 
sight is smaller than unity at all times for all three runs, 
even the supercritical one, B2-AD. This is likely a conse- 
quence of mass conservation again (cf. Sec. 14. 2p . Because 
these statistics are weighted by area over the circular re- 
gion, the supercritical regions, which are denser, occupy a 
smaller fraction of the surface area of the cloud, and there- 
fore the area-weighted average M2FR is smaller than unity 
in all cases. However, the averages for the high-column den- 
sity LOSs, shown in the right panels of Fig. [TT] are larger 
than unity for all times in run B2-AD, and for nearly 20 
Myr in run B3-AD. Instead, for run B4-AD, even this high- 
A'' average barely reaches values larger than unity, and only 
over less than 10 Myr. Finally, in all cases, the high-A*" aver- 
age M2FR decreases at late times, a phenomenon that can 
be attributed to the high-M2FR gas consumption by star 
formation. 

In general, we conclude from Sections 14.21 and 14.31 that 
the M2FR is a time-dependent function of time as a cloud is 
built up by converging WNM streams, whose average gener- 
ally increases with time, although in the high-A regions the 
average later decays due to consumption by SF. Spatially, 
the M2FR exhibits large fluctuations, whose 3-a range spans 
over one order of magnitude even in the strong-fleld cases. It 
is the high-M2FR tail of the distribution that is responsible 
for star formation. 



4.4 Buoyancy of subcritical regions 

An unexpected feature we have observed in these simula- 
tions is that the subcritical and supercritical regions do not 
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maintain their relative positions fixed tliroughout the evolu- 
tion of the simulations. Instead, the subcritical regions ex- 
hibit "buoyancy" , so that they tend to separate themselves 
from the global contraction of the clouds, even in the glob- 
ally subcritical cases that rebound at later times. This is 
most clearly seen in animations of the simulations, but can 
be seen in Fig. 1121 where we show that, as time proceeds, the 
subcritical regions become segregated from the supercritical 
ones, moving outwards, and developing cometary shapes, 
with their heads pointing outwards as well. 

This behavior can be described as a macroscopic-scale 
analogue of the very process of AD. In the latter, the neu- 
trals sink into the gravitational potential well, percolating 
through the ions, which remain attached to the magnetic 
field lines. In our case, the supercritical regions sink into the 
potential well, while the subcritical ones remain in the outer 
parts of the well, being held up by the magnetic tension. The 
process is also re miniscent of the interchang e mode of the 
Parker instability (|Hughes fc Cattaneolll987l '). 



5 SUMMARY AND CONCLUSIONS 

In this paper we have presented a study of the forma- 
tion and evolution of GMCs by the convergence of WNM 
streams, or "inflows", in the presence of magnetic fields 
and AD. As described b y many groups (see the review by 
IVazquez-Semadenill201o| . and references therein), this pro- 
cess involves the transition of the atomic gas from the warm, 
diffuse phase to the cold, dense one, allowing the fresh cold 
gas to quickly become self-gravitating and begin to contract. 

We first reviewed the general evolution of the gas and 
the M2FR expected in this t ype of systems, n oting that, 
as originally pointed out by iHartmann et al.l (|200^, the 
mass-to-flux ratio (M2FR) of the cloud is expected to in- 
crease in time, so that the cloud becomes molecular, self- 
gravitating, and magnetically supercritical at roughly the 
same time, provided that there is enough mass in the con- 
verging streams to render them supercritical. This condition 
requires that, for solar neighborhood conditions, the inflows 
extend beyond the accumulation length given by eq. ([2]). 
Flows that do not extend to such distances are expected 
to form subcritical clouds which, however, may be predomi- 
nantly atomic. This suggests that, in particular, the converg- 
ing flows induced by the spiral-arm potential wells, which 
have typical size scales ~ 1 kpc, will in general induce the 
formation of magnetically supercritical molecular clouds. On 
the other hand, converging flows induced by smaller-scale in- 
flows, such as supernova shocks, or simply turbulent random 
motions in the gas, may lead to the formation of subcritical, 
partially atomic clouds. 

We then discussed the difficulties inherent to measuring 
the M2FR, even under controlled conditions such as those 
of the simulations. Other than measuring the M2FR along 
magnetic flux tubes, which is impossible to perform obser- 
vationally, and perhaps even numerically, we considered two 
methods for measuring the M2FR: the "projection" method, 
which mimicks the observational procedure of measuring the 
ratio of column density to field strength along the line of 
sight (LOS), and which gives an upper limit to the actual 
M2FR, and the "local" method, which simply measures the 
mass and magnetic field in a grid cell in the simulation, and 



estimates the M2FR as given by eq. ((6)1 , giving a lower limit 
to the M2FR. 

We next studied the evolution of the M2FR and the 
star-forming properties of clouds formed by both sub- and 
supercritical inflows. We concluded that in our simulations, 
the effect of numerical diffusion is at a comparable level to 
that of AD. We found that the subcritical cases do undergo 
an initial phase of contraction, followed by a re-expansion, 
settling into an oscillatory regime. The supercritical case, on 
the other hand, proceeds directly to collapse, as expected. 
All cases form stars, although at greatly different rates, 
producing what appears more a continuum of star-forming 
regimes as the mean magnetic field strength is varied, rather 
than a bimodal regime of high SFR in supercritical cases and 
low SFR in subcritical ones, as was the ca se in the "stan- 
dard" model of ma gnetically regulated SF l|Shu et al.lll987l : 
iMouschoviasI [l99ll ) . In particular, the marginally subcriti- 
cal case B3-AD, through the action of diffusion, reached SF 
efficiencies (SFEs) comparable to those of the supercritical 
case B2-AD. The onset of SF is delayed by up to 15 Myr 
in the most strongly magnetized cases we studied, although 
this delay is reduced by a few to several Myr when AD is 
included (when it is not, all SF activity is due to numerical 
diffusion) . The SFEs observed in our simulations range from 
~ 35% for run B2-AD to ~ 3% for run B4-AD. However, 
since stellar feedback, which would reduce the SFE even 
further, is not included in these simulations, it is likely that 
the efficiency of run B4-AD is actually too low in comparison 
with observed values. 

We then investigated the spatial and statistical distri- 
bution of the M2FR, finding that this is a highly fluctuating 
quantity. The fragmentation of the cloud by the combined 
action of thermal, nonlinear thin-shell, and gravitational in- 
stabilities leads to the formation of clumps of high density 
and high M2FR, and of low-density, low-M2FR patches. The 
fluctuations in the M2FR we observed span between one- 
and-a-half and two orders of magnitude, the distribution 
being wider for weaker magnetizations. These results are 
qualitatively consistent with recent observational determi- 
nations suggesting that the magnetic field s trength in MCs 
is strongly fluctuating (|Crutcher et al.|[2010l ). 

Next, we discussed the evolution of the mean M2FR and 
its 3-(T range in the various simulations, finding that in gen- 
eral it initially increases as expected by the assembly process 
of the cloud, to later reach a roughly stationary regime. The 
M2FR of the high-column-density LOSs, on the other hand, 
tends to decrease at later times, due to the consumption of 
gas in this regime by SF. 

Finally, we reported the occurrence of an unexpected 
effect: the buoyancy of the low-M2FR regions with respect 
to the high-M2FR ones, through a process that appears as 
the macroscopic-scale analogue of AD: the high-M2FR re- 
gions sink deeper into the potential well of the cloud, while 
the low-M2FR ones remain supported by the field in the 
outer parts of the cloud, so that the clouds evolve towards a 
segregated state with low M2FR in their periphery and high- 
M2FR towards their centre, even on scales much larger, and 
densities much lower than, those directly affected by AD. 
The process also bears resemblance with the intercheange 
mode of the Parker instability. 
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Figure 1. Cross section through the simulation centre on the {x,y) plane, showing the density field on the plane and the component of 
the magnetic field on this plane, for the B3-MH {left panel) and for the B4-MH {right panel) runs at t = 20 Myr. The initial flow collision 
occurred along the x direction, which is also the direction of the mean magnetic field. 
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Figure 2. Snapshots of run B2-AD showing the column density integrated along the central 20 parsecs of the simulation along the 
a;-direction (perpendicular to the colliding inflows), at times t = 10 Myr {left panel), t = 20 Myr {middle panel), and i = 30 Myr {right 
panel) . 
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Figure 3. Snapshots of run B3-AD showing the column density integrated along the central 20 parsecs of the simulation along the 
a;-direction (perpendicular to the colliding inflows), at times t = 9.6 Myr {left panel), t = 26 Myr {middle panel), and t = 29.7 Myr {right 
panel). The simulation contracts until t = 26 Myr, after which it begins to rebound. 
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Figure 4. Snapshots of run B4-AD showing the column density integrated along the central 20 parsecs of the simulation along the 
x-direction (perpendicular to the colliding inflows), at times t = 9.6 Myr {top left panel) t = 20.5 Myr {top right panel), i = 34 Myr 
{bottom left panel), and t = 48 Myr {bottom right panel). The simulation is seen to be undergoing global oscillations, with star formation 
(indicated by the dots) scarcely occurring only during the times of maximum compression. 
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Figure 5. Evolution of the total dense gas mass (n > 100 cm~^) and total sink mass (top le]t panel), the time derivative of the total 
sink mass Mginkss (approximately giving the SFR; top right panel), the total number of sink particles {bottom left panel), and the SFE, 
defined by eq. (0 {bottom right panel) in all simulations. The density range for the dense gas includes gas that would classify as "atomic" 
as well as "molecular" . The graphs of Msi^j^s have been smoothed by averaging ovser 16 neighboring data points at each plotted value. 
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Figure 6. Snapshots of the mass-to-flux ratio /i, normahzcd to the critical value, for runs B2-AD {top left panel), B3-AD {top middle 
panel), B4-AD {top right panel), B3-MH {bottom middle panel) and B4-MH {bottom right panel) at time t = 20 Myr. The M2FR is 
calculated as N/Bi^os> indicated by eq. lO, integrated over a 20-pc path centered at the midplane of the simulation along the direction 
of the inflows. This path completely encloses the cloud along the x-direction in all simulations. Note that there exist vast numbers of 
supercritical filaments intermixed with subcritical patches. The latter, however, occupy most of the volume. The dots indicate the 
positions of sink particles. Note also that the density structure is similar at large scales but different in the small-scale detail between 
the cases with and without AD. 
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Figure 7. Left: Comparison of the M2FR histograms at 4 = 20 Myr, for cases with {dotted lines) and without {solid lines) AD, for runs 
B3 {left panel), and B4 {right panel). The M2FR is estimated through the "projection method", as indicated by eq. JSjl over a circular 
region centered at the point {y, z) = (0, 0), and of radius 20 pc. 
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Figure 8. Top left: Histograms of the M2FR, estimated through the projection method, for runs B2-AD {solid line), B3-AD {dotted 
line), and B4-AD {dashed line), at t = 20 Myr in all cases. Top right: Cumulative probability distributions for the same three runs, with 
the same line coding. Bottom left Density-weighted probability distribution, giving the fraction of mass below the indicated value of //. 
Bottom right Same as the bottom left panel but for high-column density LOSs {N > 10^^ cm~^) only. 
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Figure 9. Left: Histograms of tlie M2FR at t = 20 Myr, estimated tlirougli the "local method", using eq. for each cell of the middle 
plane, x = 0, of the three simulations, B2-AD {solid line), B3-AD {dotted line) and B4-AD {dashed line). These histograms can be 
compared to those using the projection method, shown in the left panel of Fig. |8] 
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Figure 10. SFE versus the supercritical mass fraction in runs B2-AD, B3-AD and B4-AD, with the SFE read off from the bottom right 
panel of Fig. \5\ and the supercritical mass fraction read off from the bottom left panel of Fig. [S] The plotted value of the SFE is the 
mean between the extremes taken by the SFE over the time interval after which the initial rapid growth has ended (18 < t ^ 26 Myr 
for B2-AD, 25 < i ^ 36 Myr for B3-AD, and 30 < t ^ 58 Myr for B4-AD), and the error bars denote these extremes. The dotted line 
indicates a least squares fit to the data points, given by eq. 
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Figure 11. Evolution of the M2FR fi, normalized to the critical value, for runs B2-AD (fop), B3-AD (middle) and B4-Ad {bottom), 
computed using the projection method over the same circular region as in Fig.|8] The solid lines show the mean value of /i//.tcrit, and the 
dotted lines delimit the 3cr range of fi, where the mean and the standard deviation are calculated for log fi. The left panels show these 
quantities computed for all lines of sight parallel to the axis of a cylindrical volume of length and diameter both equal to 20 pc, centered 
in the centre of the numerical box. The right panels show the same quantities computed only for lines of sight having column densities 
N > 10^^ cm~'^. In all cases, the lines of sight extend over the interval —10 < x < 10 pc. 



Molecular Cloud Evolution IV: Magnetic Fields 21 




Figure 12. Snapshots of the mass-to-flux ratio /i, normalized to the critical value, for run B4-AD at times t = 20 (left panel), 22 {middle 
panel), and 24 {right panel) Myr. Note the evolution of the subcritical regions, shown in green, which develop cometary shapes, pointing 
outwards from the centre of the image. 



